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Abstract 

In this work we study the transition from normal to anomalous diffusion of 
Brownian particles on disordered potentials. The potential model consists of a 
series of “potential hills” (defined on unit cell of constant length) whose heights 
are chosen randomly from a given distribution. We calculate the exact ex¬ 
pression for the diffusion coefficient in the case of uncorrelated potentials for 
arbitrary distributions. We particularly show that when the potential heights 
have a Gaussian distribution (with zero mean and a finite variance) the diffu¬ 
sion of the particles is always normal. In contrast when the distribution of the 
potential heights are exponentially distributed we show that the diffusion coeffi¬ 
cient vanishes when system is placed below a critical temperature. We calculate 
analytically the diffusion exponent for the anomalous (subdiffusive) phase by 
using the so-called “random trap model”. We test our predictions by means 
of Langevin simulations obtaining good agreement within the accuracy of our 
numerical calculations. 

Keywords: 


1. Introduction 

Transport of classical overdamped particles in random media has been in¬ 
tensively studied due to its relevance in several physical systems 0]. Due to 
its importance, several kind of models of transport in disordered media has 
been introduced in order to understand the underlying mechanisms originat¬ 
ing the observed phenomenology 00- Depending on the model type, the 
system can exhibit normal or anomalous transport accompanied of normal or 
anomalous diffusion. For example, recent works have revealed that a system 
of overdamped particles on a random potential the unbiased diffusion is always 
normal independently of the correlations 0, H]. In contrast, in models with 
uncorrelated Gaussian random force field exhibit anomalous diffusion when the 
system is unbiased (i.e., when the system is not subjected to a constant driving 
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force and consequently it does not exhibit net transport) ®S®. Particularly, 
some attention has been addressed to models which exhibit a transition from 
normal to anomalous diffusion as a function of the temperature. For example, a 
simplified model exhibiting this feature is the so called random trap model [7]. 
In this model the phase space consist of a “chain” of traps, indexed by i £ Z 
with depth AEi. The time that a particle spent inside a trap has a distribution 
whose first moment diverges at a critical temperature T c > 0. This is a result 
that follows from the statistical properties of the energy depths AEi , whose 
distribution is exponential. On the other hand, it has also proposed a model of 
Brownian particles with continuous disorder (a Gaussian squared potential) in 
which a dynamical transition also occurs at a finite temperature [8}. 

In this work we study the particle-polymer model, a model which is char¬ 
acterized by having continuous dynamics but discrete disorder @,0. We will 
show that this type of models, which can be considered in the middle between 
the random trap models (which are discrete) and models with gaussian (con¬ 
tinuous) disorder, are able to exhibit a transition from anomalous to normal 
diffusion as a function of the temperature. The particle-polymer model consists 
of an ensemble of Brownian particles interacting with a random polymer. The 
polymer is built up by concatenating at random (by means of some stochastic 
process) monomers of constant length which can be taken from a finite or infinite 
set of monomer types. The interaction of a given particle with the monomers 
defines the potential profile that the particle sees on a specific unit cell. When 
the interaction is short-ranged we can think of this model as an ensemble of par¬ 
ticles moving on a series of tracks, each track having a fixed random potential 
profile. Based on the theory developed in Ref. 0 , we will show that when the 
potential profile has barrier heights with normal distribution, the diffusion is al¬ 
ways normal. In contrast, when the distribution of the potential barrier height 
is exponential, we will show that there is a finite temperature T c at which the 
system transits from normal to anomalous diffusion. We use the random trap 
model to approximate the behavior of the system in order to obtain an an an¬ 
alytical expression for the diffusion exponent in the anomalous (subdiffusive) 
phase. 

In consequence this work is organized as follows. In section [2] we introduce 
the model that we study and review some results already known for this system. 
In Section[3]we calculate the diffusion coefficient for the systems by means of the 
Einstein relation which is known to be valid for our model. We also calculate 
explicitly the diffusion coefficient for the cases in which the potential height has 
a Gaussian and exponential distributions respectively we show that in the latter 
case the diffusion coefficient is zero below a critical temperature. In Section [4] 
we study the case in which the potential heights are exponentially distributed 
and the temperature is below the critical one. Since in this case the diffusion 
coefficient is zero the system exhibit anomalous subdiffusion. We implement 
what we call the “random trap” approximation. This approximation is based 
in identifying the waiting time distribution for the random trap model with 
the distribution of the mean first passage time induced by the disorder. This 
method allows us to obtain a good approximation for the diffusion exponent for 
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the subdiffusive anomalous phase. We compare all the predictions made for this 
system against numerical simulations. Finally in Section [5] we give the main 
conclusions of our work. 


2. The model 

The particle-polymer model for transport in disordered media has been put 
forward in Ref. [9]. This model consists of an ensemble of non-interacting over¬ 
damped particles sliding over a ID substrate. The substrate is model as a series 
of unit cells of constant length, called monomers if the substrate is interpreted as 
a polymer, with which the particles are interacting. This model mimics in some 
sense, the motion of a given protein along a biological “disordered” polymer 
such as the DNA. The motion of a given particle along the referred substrate is 
ruled the stochastic differential equation, 

ldX t = (f(X t ) + F)dt + Qo dW t . (1) 

Here X t represents the particle position at a time t, f(x) is minus the gradient 
of the potential V(x), which is the potential induced by the particle-polymer 
interaction and W t is a standard Wiener process. The constants q$, F and 7 
are the noise intensity, the external driving force, and the friction coefficient 
respectively. According to the fluctuation-dissipation theorem Qq = 2y/3 —1 , 
where /3, as usual, stands for the inverse temperature times the Boltzmann 
constant, /3 = l/fc^T. 

Of course, the interaction of the particle with the polymer specifies the po¬ 
tential profile seen by the particle when it is located at a specific monomer. 
If the particle-polymer interaction is short-ranged we can assume that the po¬ 
tential profile seen by the particle at a given monomer is defined exclusively 
by the monomer type at which the particle is located. Let us denote by 
a :=(... , a_i, ao, ai,...) a specific realization of the random polymer, where a, 
is the monomer type at the jth unit cell. The variable aj can be interpreted 
as a random variable defining the interaction with the particle or defining the 
properties of the potential profile if the interaction is short-ranged. Thus, in the 
latter case, we can think of the potential V(x) as a function of both, the particle 
position x and the random variable a n , where n is index identifying the unit 
cell at which the particle is located. If the particle is located at the nth unit 
cell, then its position can be written as x = nL + y , where L is the length of the 
unit cell and y is the relative position of the particle on the monomer. Thus, 
the potential V ( x ) can be written as a function i/j of two variables: the relative 
position y and the random monomer type, V(x) = '4>{y 1 a n ). This model for dis¬ 
ordered media was first studied in absence of noise in Ref. [1] and the authors 
showed that the particle current and the effective diffusion coefficient can be 
known exactly for a large class of disordered potentials. Later on, in Ref. 0 it 
was shown that the particle current and diffusion coefficient can still be exactly 
in terms of quadratures. However in Ref. (lo| the author was mainly interested 
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in the case of biased transport, and, in fact, no formula for the diffusion coef¬ 
ficient in absence of driven force was done. Here we will explore the unbiased 
diffusion based on the findings of Ref. jToj. 

First of all we should remind that in Ref. [10] the author gave an exact 
expression for the first passage time (FPT) on a unit cell, averaged with respect 
to the noise, given by 


OO 


Ti(a)= 7 /3^ 

777 .= 1 

e- m ^ i g + (a) g _[a- m (a)]+ 7 /3/o(a). 

(2) 

In the above equation a stands for the realization of the polymer and q± and Iq 
are integrals defined as 

<?±(a) := 

/ dxexp ( zb /3[^(x, a) — xF ]), 

Jo 

(3) 

4(a) = 

/ Q~(x, a)B + (x, a)dx, 

Jo 

(4) 


where the functions B± : R x y —* R and Q± : R x _4 Z —> R were defined as, 


(x, ci) 

:= / dyexp(±P[ijj(y,a)-yF]), 

do 

(5) 

B±(x, a) 

:= exp ( ±/3[^>(a;,a) — xF]). 

(6) 


The above expressions are valid for a large verity of disordered potentials, 
since these expressions were derived under some mild general assumptions. 
Specifically, the two restrictions imposed were i) the ergodicity of the random 
polymer when the sequence of monomer is seen as a discrete-time stochastic 
process and ii ) the validity of the central limit theorem in its general form [ill 
ll2l . 1131 1 • These two properties can be exploited in order to consider more general 
random potentials, particularly those that can give rise anomalous diffusion in 
absence of external driving forces. First, let us remind that the particle current 
can be written in terms of the FPT, averaged with respect to the noise and the 
polymer ensemble, as 0, 

Jes(F) = -i-, (7) 

where we used the notation (-) p for the expected value with respect to the 
polymer ensemble. Particularly, when the Brownian particle has a short-range 
interaction with the polymer, we can assume that the potential only depends on 
one “coordinate” of the polymer (as we explained above) and that the monomers 
along the polymer does not have correlations, i.e., the random polymer can be 
considered as a result of a realization of an infinite set of independent and 
identically distributed (i.i.d.) random variables. These hypotheses allow write 
down a more explicit expression for (Ti) p as follows, 

—/3F l 

(Ti(a)) p = 7/3 1 _ (g+(a))p(g-(«))p + bff (4(«)) p ■ (8) 
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Notice that in the last relation we have no longer made explicit the dependence 
of T\ with the whole polymer, since the potential depends on only one monomer, 
and, due to translational invariance, the position of such a monomer does no 
longer matters. Thus, the polymer average will only depend on the random 
variable a and its distribution. 


3. The diffusion coefficient 


It is known that in the class of disordered systems we are considering the 
Einstein relation remains valid [l|. This allows us to obtain an expression for 
the diffusion coefficient in absence of driving force. If we denote by _D e ff the 
effective diffusion coefficient for the system in absence of external driving force 
(i.e., for F = 0), the Einstein relation states that 


D e ff = s lim 
/3 F^o 


Jeff(F) 

F 


To obtain an expression for the diffusivity in this case, first notice that T\ 
diverges as F goes to zero. It is easy to see that the asymptotic behavior is 
giving by, 

(Ti(a )) p « J^Z^Zir-P) + 7 /? </ 0 (a)) p , 
where Z(0) is defined as, 


m 



From the above expression we can see that the diffusion coefficient can be written 
straightforwardly as follows, 


DeS 1 pz(pz(-py 

which is the effective diffusion coefficient at zero driving force, F = 0. Below we 
will use this expression to exactly calculate the diffusion coefficient in a specific 
potential model. 


3.1. The potential model 

The model potential that we will use consists of a piecewise linear function. 
This makes the force field to be piecewise constant. Let x be the particle position 
and let n be the cell at which this particle is located. It is clear that x can be 
written x = y + nL, where y is the relative position of the particle on the nth 
unit cell. Thus, the potential that the particle feels at x is defined as, 


V(x) = ip{y,a n ) 


2a n y/L if 0 < y < L/2 
2a n (L — y)/L if L/2<y<L. 


(9) 
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where a:= (..., a_i, ao, ai,.. .) represents a realization of the random polymer. 
Here a n stands for the potential height and this quantity specifies a polymer 
type (see Fig. |T]for a schematic representation of the potential model). Thus 
the set {..., a_i, ao, ai,...} can be considered as a realization of an infinite 
set of i.i.d. random variables. The reason to take these random variables as 
independent is because the polymer is assumed to be uncorrelated. On the 
other hand, the reason to assume that the random variables are identically 
distributed is just to fulfill the condition of translational invariance. 



Figure 1: Schematic representation of the potential model, (a) The potential profile on the 
Oth unit cell, (b) A realization of the random potential with a few unit cells. The monomer 
type is labeled by a n , with n £ Z and this label is in turn a random variable giving the height 
of the potential hill. There are as many monomer types as the number of values that can take 
any random variable a n . 

Let a be a random variable with the same distribution of any a n for n £ Z. 
If the random variable a is assumed to be discrete taking values on a finite set 
W ■■ o < j < TV}, then, the number of monomer types is finite. However, the 
random variable a is not restricted to take values in a finite set nor a discrete 
one. The random variable can even take values in a continuous set and the 
hypotheses of polymer ergodicity and translational invariance is not violated. 
In this case the number of monomer types will be uncountably infinite. It 
is particularly interesting to investigate the case in which the set of possible 
potential heights is R or R + , i.e., when the height distribution have unbounded 
support. In this case, we are in a situation in which the potential hills have 
unbounded height and the particle could stay long times in a valley along the 
polymer. In this situations we would expect to have a competition between 
the thermal activation, taking the particles out of the valley, and the height of 
potential barriers, slowing down the dynamics of the particles. Below we show 
that to slow down the diffusivity is such a way that the diffusion coefficient 
be zero it is not enough that the height distribution have unbounded support. 
Some additional features on the distribution are needed in order to have a zero 
diffusion coefficient. In particular we show that when the height distribution 
is gaussian, with zero mean and finite variance, the diffusion coefficient never 
vanishes. On the other hand, when the height distribution is exponential, the 
diffusion coefficient can be zero below a critical temperature, i.e., the system 
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freezes” at a finite (positive) temperature. 


3.2. Gaussian heights 

Let us analyze the case in which the height distribution is Gaussian with zero 
mean and variance a 2 . The variance a 2 is an intrinsic parameter characterizing 
the disorder. Along this line we can say that the probability density function 
p(a) of the height of the potential hills along the polymer is given by, 

PgO) = exp (-« 2 /2 o- 2 ) • 

V 2tt(j z 

Note that the distribution allows negative values for the potential “height” which 
mean that a typical realization of the potential contains not only “potential hills” 
but also “potential wells”. 

In order to calculate the exact diffusion coefficient we should first evaluate 
Z((3 ), which can be written as 



dxex p f3ip(x,a)]p g (a). 


First notice that the potential ip(y,a n ) can be factorized as the random height 
a n times a dimensionless function tp which is piecewise linear, i.e., ip(y,a n ) = 
a n ip(y), where 


<p(y) = 


2y/L if 0 < y < L/2 
2 (L - y)/L if L/2<y< L. 


( 10 ) 


Then it is easy to see that the integral Z(/3) can be written as , 

pL /»oo 

Z(/3) = dx daexp[/3aip(x)] exp (— a 2 /2a 2 ), 

J0 J —oo V27TCT 2 

where the integral with respect to the variable a can be easily evaluated to give, 

Z(/8) = f dxex p \p 2 a 2 p 2 [x)!2\ . 

Jo 

Now, the above result allows us to calculate Z{fJ) for our model explicitly in 
terms of quadratures. After a few manipulations it is possible to see that, 

Z(J3) = 

po- 


where we have defined G as the function, 

(■po¬ 
lo 


j-po 

G{f3a) = / exp (x 2 /2)dx. 
Jo 
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( 11 ) 


The result for Z(P) allows us to write the diffusion coefficient as 

D = ^ 

eff r yG 2 (pa) ’ 

Notice that when the temperature is large, i.e., when pa -C 1, we can ap¬ 
proximate the function G{pa) as, 

G(/3a) « Pae 13 a ! 2 for pa <C 1. 


Then, the diffusion coefficient take the form, 

—/3 2 ct 2 

D e s ~ -—jj— for Pa < 1, (12) 

7 P 

which 

disord □ □□ 



Figure 2: Effective diffusion coefficient for the random potential with Gaussian height dis¬ 
tribution. We display the diffusion coefficient given by the exact expression of Eq. m (solid 
line) and by the numerical experiment (filled circles). For the numerical simulation we gen¬ 
erate five polymer of 2 x 10 6 monomers. On each monomer we assumed that the potential 
profile has the form of a “hill”, which is given by Eq. J10D . The potential heights {a n } were 
defined by generating a series of 2 x 10 6 normally distributed random number with zero mean 
and variance a 2 = 1. Then we placed 4000 particles along the polymer which were subjected 
to the Langevin dynamics simulations. The obtained trajectories were used to estimate the 
mean square displacement, which in turns gave us an estimation for the effective diffusion 
coefficient. 

In Fig. [2] we display the effective diffusion coefficient analytically calculated 
using Eq. CD- In the same figure we also display the effective diffusion coef¬ 
ficient obtained from the numerical simulation of an ensemble of overdamped 
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particles in a Gaussian random potential. The random potential was obtained 
by generating five series of 2 x 10 6 Gaussian random numbers with zero mean 
and variance one (a 2 = 1). Then, we used these random numbers as the heights 
of the potential hill comprising the polymer. Thus every polymer has 2 x 10 6 
monomers. Next, we placed 4 000 particles equally spaced on every polymer and 
such particles were subjected to the Langevin dynamics during a total time of 
10 5 arb. units. We fixed the friction constant 7 equal to one. Then we used the 
resulting trajectories to estimate the mean square displacement, from which we 
obtained the corresponding estimation of the diffusivity. As we can see in Fig. [ 2 ] 
w have a good agreement (within the accuracy of our simulations) between the 
exact curves an the numerical experiment. 


3.3. Exponential heights 

We now consider the case in which the height distribution is exponential. 
Notice that for the exponential distribution there is no negative heights, which 
means that there is no potential “wells” but only potential “hills”. We could 
think that, since there is no potential wells the probability that a given particle 
gets stuck is, in this case, lower than in the case of the Gaussian height distribu¬ 
tion. This is because in the present situation the particle has to overcome only 
the potential barrier due to the hills, contrary to the Gaussian case in which 
the particle has to overcome the potential barrier when it is located at a well 
(negative height) plus the potential barrier due to the presence of a hill (positive 
height). However we will show that, contrary to our intuitive arguments, the 
diffusivity in the exponential case is lower than in the Gaussian case. This slow¬ 
ing down of the diffusivity is still much more dramatic, since at sufficiently low 
temperatures the diffusion coefficient vanishes. This means that there occur a 
dynamical transition from normal to anomalous subdiffusion as the temperature 
decreases. Moreover, we will also show that the mean square displacement of 
the Brownian particles grows in time asymptotically as AX 2 ~ t a where a, the 
diffusion exponent, can be analytically calculated by using the approximation 
provided by the random trap model for subdiffusion. 

First let us calculate the diffusion coefficient for this model. We assumed 
that the height distribution is exponential, which means that the corresponding 
probability density function reads as, 


Pe(a) 


exp (—a/A) 
A 


As in the Gaussian case we need to calculate the function Z(f3) which in this 
case can be written as, 


Z{P) 


OO nL 

/ da dxexp [P'i/;(x,a)]p e (a) 
1 0 Jo 

pL poo 


dx / da exp [/3ay>(a;)] 


'0 Jo 

pL 


exp (—a/A) 

A ' 


dx 


1 0 1 - P^p{x) 
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Notice that the last integral diverges when pXip(x) = 1 at some x € [0, L\. 
This is of course possible for temperatures low enough. Indeed if we recall the 
definition of ip [see Eq. (HI] we observe that 0 < <p(x) < 1 for all x G [0, L\ and 
in particular the upper bound is reached at x = L/ 2. This implies that when 
the critical (inverse) temperature is given by p c = 1/A. If P > j3 c then Z(/3) 
diverges and consequently the diffusion coefficient is zero. Conversely, if /3 < /3 C 
then Z(/3) is finite and the diffusion coefficient is positive. In the latter case the 
integral can be done exactly for our potential model, giving, 

Z{P) = ~ ln(l - PX) for p < p c = 1 /A. 

The above result allows us to write down explicitly the effective diffusion coef¬ 
ficient, 


D e fi 


px 2 1 

7 |ln(l-/3A)ln(l + /?A)| 


for P < p c = 1/A. 


(13) 


To perform the simulations of the Langevin dynamics as well as to numeri¬ 
cally evaluate D e g let us introduce dimensionless quantities. First notice that p\ 
is a dimensionless parameter which allows us to define naturally a dimensionless 
temperature T as 


T := 


1 

PX' 


Similarly we should observe that A /7 has the same units as the diffusion coef¬ 
ficient. This suggest to define a dimensionless effective diffusion coefficient D e g 
as 


-D e ff := 


7-^eff 

A 


Using the above definition we can rewrite the expression for the diffusion coef¬ 
ficient in terms of dimensionless parameters, 


D e fi 


T 


i»(l-i)ln(l + i) 


for f > T c = 1. 


In Fig. [3] we display the effective diffusion coefficient above the critical tem¬ 
perature. To perform the numerical simulations we have chosen the parameters 
A = L = 1. As we saw above we have that the dimensionless critical temper¬ 
ature T c equals one. We observe that within the accuracy our simulations the 
numerically obtained diffusion coefficient agrees with the one obtained by means 
of the exact closed formula given in Eq. m- 


4. Anomalous subdiffusion: an analytical approach 

We have seen that when the heights are exponentially distributed the dif¬ 
fusion coefficient goes down to zero for all the temperatures below a critical 
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Figure 3: Effective diffusion coefficient for the random potential with exponential height 
distribution. We display the diffusion coefficient given by the exact expression of Eq. m 
in terms of dimensionless quantities (solid line) and by the numerical experiment (filled cir¬ 
cles). For the numerical simulation we generated five polymer of 2 x 10 6 monomers. On each 
monomer we assumed that the potential profile has the form of a “hill”, which is given by 
Eq. (1101 ). The potential heights {a n } were defined by generating a series of 2 x 10 6 exponen¬ 
tially distributed random number with zero mean and variance cr 2 = 1. Then we placed 4 000 
particles along the polymer which were subjected to the Langevin dynamics simulations. The 
obtained trajectories were used to estimate the mean square displacement, which in turns 
gave us an estimation for the effective diffusion coefficient. 


temperature T c = X/ks- This fact implies that the mean square displacement 
AXt grows in time lower than linear, i.e., 

AA' t 2 ~ t a (14) 

where a is what we call the diffusion exponent. The diffusion exponent charac¬ 
terizes the subdiffusive regime, indicating in some sense how low is the diffusivity 
of the Brownian particles on the random potential. Here we will implement an 
analytical treatment in order to obtain an approximate expression for the diffu¬ 
sion exponent a. First of all we should notice that the origin of the subdiffusive 
behavior lies on the statistics of the potentials. The height of the potential 
barriers are responsible of the (long) time that a particle remains in a given cell. 
The higher potential barrier the larger time the particle spent on the cell. 

Thus we need to estimate the distribution of the time that a particle takes 
to go from a given cell to the adjacent one. This time will be latter identified 
as a “waiting time” in a random trap model. To perform such an estimation we 
assume that the time that a particle need to go from a unit cell to the adjacent 
one is given by the inverse of the Kramers escape rate K esc . This quantity is 
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given by [? ] 


ft-esc — 2 ^ \/ V n (^min) V ,f (^max) GXp ( /5 [V(^max) ^(#min)] ) • 


It is important to remark that the Kramers escape rate is an approximation 
to estimate the time that a particle need to jump to the next unit cell. This 
is because the “true” waiting time (let us call T w ) is a random variable that 
depends on the two underling stochastic processes present in the system: the 
“Gaussian white noise” (i.e., the Wiener process W t ) and the random polymer 
a. We actually pretend that waiting time is nearly equal to the inverse Kramers 
t := n~ s ], escape rate, which is asymptotically equal only in the limit of zero 
noise strength. The latter is due to the fact that is in fact the average of 
T w with respect to the noise, i.e. r = (T w ) p . Thus we naturally have that t is 
a random variable that only depends on the realization of the random polymer 
a. Moreover, since the potential model we have introduced only depends on one 
coordinate of the random polymer, we have that r only depends on one random 
variable (or, in other words, r only depends on one coordinate or a). 

In this case it is easy to see that the escape rate is approximately given by, 

k ~ ^ r -P a 

esc 2 ttL 2 ’ 

which means that the waiting time can be approximated as 


r(a) 


2t tL 2 
X 


p/3a 


(15) 


A direct calculation shows that the distribution function for r can be written 


as 


fr(t) 


1 / A t \ _( ^ +1) 
2tt/3 V 27 tL 2 ) 


for t large enough. 

Then, following the assumption that the distribution of r corresponds to 
the waiting time distribution of a random trap system, we have, according to 
the theory developed in Ref. 11], that the diffusion exponent a defined through 
Eq. (Il4l) is given by, 


a = 


2 

1 + /3A' 


(16) 


In terms of the dimensionless quantities, defined above, a can be rewritten as, 


a = 



(17) 


In order to verify the approximations made to obtain the analytical expres¬ 
sion for the diffusion exponent we have performed numerical simulations. Firstly 
we made “random polymers” by producing five series of 2 x 10 6 exponentially 
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Figure 4: Diffusion exponent for the random potential with exponential height distribution. 
The figure shows the diffusion exponent obtained from the approximate analytical expression 
of Eq. (1171 ) (solid line) and from the numerical experiments (filled circles). The numerical 
experiment was done by simulating the stochastic dynamics of 20 000 particles, every particle 
on a different realization of the polymer. These trajectories were used to obtain the mean 
square displacement as a function of time. The latter is assumed to behave as t a and from a 
fit of least squares we estimated the diffusion exponent a. 


distributed random numbers. Each series represents a random polymer of 2 x 10 6 
symbols long. Then we placed 4 x 10 3 particles evenly spaced along every poly¬ 
mer. Next the particles were subjected to the Langevin dynamics, from which 
we estimated the mean square displacement as a function of time. This proce¬ 
dure yields the effective diffusion coefficient which is displayed in Fig. [I] As we 
see, the behavior of the analytically calculated diffusion exponent is consistent 
with the simulations. This result allows us to say that the random trap model 
is a good approximation to the problem of an overdamped particles in a random 
exponential potential in presence of thermal noise. 

5. Discussion and conclusions 

In this work we have studied the particle-polymer model for diffusion on 
random potentials. We considered an ensemble of overdamped particles moving 
on a one dimensional random substrate which we interpret as a polymer. The 
substrate was assumed to interact with the particles in such a way that the 
substrate can be thought as a series of tracks (unit cells) each track having a 
fixed potential profile. We assumed that the potential height on each track is a 
random variable and we explored the cases in which the potential heights have 
Gaussian distribution and exponential distribution. We calculated the exact 
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diffusion coefficient in both cases and we found that when the height distribu¬ 
tion is Gaussian the diffusion is normal at all temperatures. In contrast, when 
the potential heights are exponentially distributed we found that the diffusion 
coefficient vanishes for low enough temperatures. For the potential model with 
exponential height distribution we were able to exactly determine the critical 
temperature T c . For temperatures below T c we assumed that the system can be 
seen as a random trap system and we estimated the waiting time distribution 
by using the Kramers escape rate. Within this approximation we were able to 
obtain the diffusion exponent as function of the temperature in a closed way. We 
performed simulations of the corresponding Langevin dynamics which confirmed 
our theoretical results within the accuracy of our numerical experiments. 
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